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When topological insulators possess rotational symmetry their spin lifetime is tied to the scattering 
time. We show that in anisotropic topological insulators this tie can be broken and the spin lifetime 
can be very large. Two different mechanisms can obtain spin conduction over long distances. The 
first is tuning the Hamiltonian to conserve a spin operator cos cj) a x + sin cj) a y , while the second is 
tuning the Fermi energy to be near a local extremum of the energy dispersion. Both mechanisms 
can produce persistent spin helices. We report spin lifetimes and spin diffusion equations. 

PACS numbers: 73.43.-f, 72.25.-b, 72.10.-d, 85.75.-d 



I. INTRODUCTION 

Topological insulator exhibit a gap in the spectrum 
of bulk states, and bridging that gap is a band of surface 
states; if the Fermi energy is within the gap then elec- 
trons flow only along the surface and not in the bulk. At 
small enough momenta the surface band has the shape 
of two cones joined at their ends. One Dirac cone de- 
scribes electrons with positive energies, and the other de- 
scribes negative energies. An electron's spin is locked to 
its momentum, so backscattering is suppressed. All these 
properties are consequences of time-reversal (T) symme- 
try, and are robust against small perturbations which are 
T-symmetric, such as non- magnetic impurities. 

Recently much attention has been given to creating 
topologically protected qubits at the interface between a 
3-D topological insulator (TI) and a conventional super- 
conductor, allowing robust quantum arithmetic^. The TI 
spin-momentum locking attracts attention to spintron- 
ics; recent works have shown that circularly polarized 
light could induce spin currents^ and topological phase 
transitions 7 . 

This paper's main focus is on obtaining good spin con- 
ductors suitable for spintronics. Disordered TIs are un- 
usually poor spin conductors. Because electronic spin is 
tied to momentum, each scattering event randomizes the 
spin, as is typical of Elliot- Yafet spin relaxation^. On the 
Dirac cone the spin lifetime r s is tightly coupled 9 to the 
scattering time r by r s /r = 2. Ordinary semiconductors 
have much longer spin lifetimes, because spin is conserved 
during scattering and is randomized only by precession 
between scattering events. In this paper we will describe 
how to tune a TI for very long spin lifetimes (r s ^> r), 
allowing spin to conduct over long distances. 

We will show that the two spin profiles shown in Figure 
[I] conduct in a properly tuned TI. The first profile lies in 
the surface plane and its angle <p in that plane is constant. 
The spin component S<j> aligned with <j> conducts: its in- 
tegral J dx 5^, = cos(j)S x +sm(j)Sy is conserved (r s = 



FIG. 1: Two spin profiles that conduct in properly tuned TIs. 
The light blue and dark blue represent the conducting surface 
and the insulating bulk of a 3-D topological insulator. Line 
a shows a spatially uniform spin density. Line b shows a 
persistent spin helix which oscillates along the x axis and is 
spatially uniform along the y axis. The PSH is a standing 
wave; it does not precess. 



oo ) and its long- wavelength variations equilibrate diffu- 
sively, similarly to heat diffusion. Spin at right angles to 
the unit vector is filtered out very quickly, relaxing 
with lifetime r s = r. Figure [TJd shows the second spin 
profile which conducts in a properly tuned TI: a standing 
spin wave which repeats at intervals of tt/\Q\. Because 
it rotates its spin orientation and does not decay, it is 
called a persistent spin helix (PSH) 10 . Associated with 
the PSH is a conserved quantity J dxSpsH, Spsh — 
cos(2Q • x)(— sin <j)S x + cos 4>S y ) + sin(2Q • x)S z . If 4> = 
the PSH rotates in the y — z plane, while if <fi = ir/2 it 
rotates in the x — z plane. The spin orientation (/>, the 
wave- vector 2Q, and the diffusion constant all depend on 
the details of the TI Hamiltonian. 

Tuning for spin conduction is not possible unless the 
surface band is more complex than a simple Dirac cone. 
Anisotropy, i.e. violation of rotational invariance^, is 
key to long distance spin conduction. This is the origin 
of the spin orientation angle (j) which parameterizes both 
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FIG. 2: Anisotropy in TIs. (a - isotropic) the Rashba Hamil- 
tonian Hr — hvp{k y a x — k x a y ), (b - stretched) linear terms 
from a k • p model of f3 — Ag^Te, (c - infinitely stretched) a 
tuned model Hi = Hvf k x (cos <j> a x + sin <j)(T y ) which conserves 
the spin operator cos0 o x ~~\~ sin cf) a y . The colored surfaces 
show the energy dispersion E(k), and the lines inscribed on 
them show Fermi surfaces - contours of constant energy. The 
arrows show the spin orientation, which is locked to the mo- 
mentum. 



spin profiles. Anisotropy can be realized in TIs eith er by 
choosing reduced-symmetry materials like (3 — Ag2T d 11 ^ 
(see Figure J2Jd) or by cutting a high-symmetry TI in a 
way that reduces the surface's symmetr ^ 13 * 14 *. Figure [2] 
illustrates distortion of the Dirac cone as rotational sym- 
metry is progressively broken. We use simple linear mod- 
els that are appropriate near the Dirac point. Figure [2^i 
shows the dispersion of a rotationally symmetric system. 
Figure [2}d shows the Dirac cone stretching along one 
axis as rotational symmetry is broken. Finally Figure 
[2J: extrapolates the stretching to its extreme: the energy 
dispersion depends only on k x not k yi and therefore an 
operator cos (j)<J x + sin <j) <j y is necessarily conserved. 

There are two ways to tune for spin conduction. The 
first is to tune the Hamiltonian for conservation of a 
spin operator cos (j)a x + sin 0cr y , causing to conduct. 
This has been achieved in GaAs quantum wells by tuning 
both the well width and the well depth to obtai n partial 
cancellation of the Rashba and Dresselhaus term d 10 * 15 * 16 *. 
These tuned quantum wells are modeled by a Hamilto- 
nian that includes both a spin-conserving quadratic term 

h 2m an d a small linear spin-orbit term which conserves 
cos (pa x + sin 0cr y , as illustrated in Figure 13b. They man- 
ifest a PSH-induced strong enhancement 01 the spin life- 
time. The TI model of Figure [2^ also conducts both S<j> 
and PSH's with wave- vector ±2Q = ±2^x. 

We will show that there is a second way to obtain 
spin conduction if the surface band has local extrema 
at k = ±Q ^ 0, in which case tuning the Fermi energy 
near the energy Eq of the extrema will produce very 
long spin lifetimes. Figure [3] shows that even a small 
anisotropy will produce the required local minima. The 
model shown in Figure [3^i includes a quadratic term and 
a spin-orbit term, both of which are rotationally sym- 
metric. It shows four Fermi surfaces - contours of con- 
stant energy. Two Fermi surfaces are at Ep = —0.163, 




FIG. 3: Fermi surfaces showing that anisotropy causes local 
,2 

extrema. H = ^— + v x k x a y + v y k y a x , m — 0.5, v y = 1. The 
solid blue lines are Fermi surfaces at energy Ep = —0.163, 
while the dashed red lines are at Ep — —0.24. All quantities 
are unitless. (a) Isotropic v x = 1, (b) anisotropic v x = 0.8, 
(c) spin conserving v x = 0. (a) and (b) are typical of untuned 
spin-orbit couplings in quantum wells, while (c) is the same 
as Bernevig et al's model 10 of a tuned quantum well. The 
dashed Fermi surfaces in (b) and (c) lie near local minima. 
Model (b) conducts spin S x when Ep ~ Eq, while (c) always 
conducts S x . The arrows point out the Fermi surface nesting 
symmetry which causes PSH's in both models. 



two are at Ep = —0.24, and all four are perfect circles. 
The minimum of the energy dispersion is also a circle 
located midway between the dotted Ep = —0.24 Fermi 
surfaces. In Figure [3}d the rotational symmetry is only 
slightly broken - the spin-orbit term's strength is only 
20% smaller along the x axis than it is along the y axis. 
Nonetheless the Fermi surfaces have already divided and 
wrapped themselves around the dispersion's local min- 
ima which are now two discrete points located inside the 
dotted E F = —0.24 Fermi surfaces. Model [3}d does not 
conserve any spin operator, but we will show that it con- 
ducts spin when Ep is adjusted so that the Fermi surfaces 
lie close to the the local minima. Lastly Figure (3^ re- 
duces the spin-orbit term along the x axis to zero, and 
exhibits spin conduction at any value of Ep. In both 
Figures [3}d and ^ spin conduction is associated with 
there being two disconnected Fermi surfaces centered on 
the local extrema. Spin conduction in TI's is possible 
only when anisotropy changes the global structure of the 
Fermi surface. 

In order to understand this global physics, we calculate 
charge and spin conduction for a very general class of T- 
conserving spin-orbit Hamiltonians: H gen = a x (k)cr x + 
a y (k)a y + ai(k), where a x: a y are odd in k and ai is even. 
This describes charge moving in the x — y plane, on one 
surface of a 3-D TI . We study the single-particle density 
matrix p = ipty, which is a 2 x 2 matrix in spin space. 
We write it as a 4— vector p = [N, S x , S y , S z ] con- 
taining the charge density N = Tr(p) and spin densities 
S x ,S y ,S z = Tr(p<Ji)/2. We will first analyze f?s struc- 
ture, and later calculate its diffusion induced by disorder. 
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II. SPIN PROFILES 



The spatially uniform and PSH profiles shown in Fig- 
ure Q] manifest themselves in the Fourier-transformed 
density matrix p(q) oc j dkip^ (q/2 + k)ip(—q/2 + k) as 
strong peaks at q = and at <f = ±2Q. These peaks are 
directly linked to the TI Fermi surface: the state vector 
ip is populated only by states from the Fermi surface, 
and therefore p(q) is peaked when q maximizes the inter- 
section of the Fermi surface with a copy of itself shifted 
by q. (The dominance of the Fermi surface is assured 
if Epr/h ^> 1, the temperature is small, and there are 
no interactions.) The spatially uniform peak p(q = 0) 
realizes this maximization trivially. 

Figures ^jp and |3j3 illustrate a special nesting sym- 
metry which produces PSH peaks at q = ±2Q. They 
show pairs of Fermi surfaces centered at ±Q which pos- 
sess inversion symmetry (Q — » —Q) because of the TPs T 
symmetry. The Fermi surfaces also possess nesting sym- 
metry, which means that a shift of q = ±2Q moves one 
Fermi surface on top of the other. This nesting symme- 
try produces peaks in p(q) at q = ±2Q; it is responsi- 
ble for PSH's. The nesting symmetry can be written as 
Q + k — > — Q + k, where Q + k lies on one Fermi sur- 
face and — Q + k lies on the other one. T plus nesting 
implies that the Fermi surface near Q possesses inversion 
symmetry around Q: Q + k — >• Q — k. 

The TI surface has only one conduction band and only 
one valence band, with spin quantum number equal to 
s = +1,-1 respectively. We assume that the Fermi 
surface lies only in the conduction band; s = +1. As 
a result the spatially uniform S z spin density is iden- 
tically zero: S z (q = 0)'s contributing terms are of the 
form S z oc Tr(pk(T z )/2, p^ — ip^ (k)ip(k), which is zero 
for all ip in the conduction band. Similarly the charge 
N(q — 2Q) and spin S<f>(q = 2Q) components of the 
PSH are zero, because p(2Q) oc ip*(Q)ip(-Q) and T 
symmetry ensures that N(2Q) = Tr(p(2Q)) and S<j> = 
Tr(p(2Q) (cos (j)cr x + sm^a y )/2) are zero when ip is in 
the conduction band. 



III. DIFFUSIVE CONDUCTION 



We now consider adding a non-magnetic "white noise" 
disorder potential V to the general Hamiltonian H gen , 

^ u(f), (u(r)u(r)) = riiU^5(r — r), 



where V 
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and riiUQ gives the disorder concentration and strength. 
When disorder is present the density matrix evolves dif- 
fusively at time scales larger than the elastic scattering 
time. Its evolution is controlled by the partial differential 
equation = 0? where the 4x4 matrix is called 

the diffuson. The diffuson's matrix structure couples the 
charge and spin densities to each other. We derive the 
diffuson using standard methods from the diagrammatic 
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FIG. 4: Diffusion diagram with two joint scatterings. A single 
joint scattering event Uj corresponds to a pair of Green's 
functions G A and G R connected by a dashed line. The dotted 
legs correspond to the external propagators. 



technique for disordered systems^K^, couched in the 
notation of References I9|201 Within this diagrammatic 
technique the conductivity is determined by the disorder- 
averaged two-particle correlation function, which is con- 
trolled by ladder diagrams at leading order in (Eft/H) -1 . 
The diffuson is composed of an infinite series of ladder 
diagrams like that seen in Figure [4j These diagrams 
describe sequences of events in which ip and ip^ move to- 
gether, scattering in unison. A single joint scattering 
event is described by the operator Z^-, and the diffu- 
son sums diagrams with any number of joint scatterings; 
Vij(q,u)) = ^2^L (Iij) n = (1 - A?') -1 - Tne j° int scatter- 
ing operator Iij is pictured in Figure [4] and is given by 
the integral 



Iij = ^J dkTr(G A (k-q/2,E F ) 
Oi G R (k + q/2, E F + Hu) aj ) 



(1) 



G A and G R are the disorder-averaged single-particle 
Green's functions which express uncorrelated movements 
of ip and ip\ while <f is the diffuson momentum. The trace 
is taken over the spin indices of G A , G R , (Ji^ and Oj, which 
are all 2 x 2 matrices in spin space. 

The zero-frequency component of the diffuson T>ij (q = 
0,cj = 0) is equal to r 8 /r, where r 3 is the tensor that 
governs relaxation of spatially uniform spin profiles. As- 
suming as before that the Fermi surface is dominant 
(Er/h ^> 1) and contains only the conduction band, we 
find: 



t s /t 





(1 - (cos2(9) F ) -(sin 20) F 
-(sin20) F (l + (cos20) F ) 



n -1 



(2) 



This result is fully general for all H gen . The angle 0(k) 
gives the relative strength of the a x cr x and a y <j y terms, 
and is defined by tan0 = a y /a x . The average ( )p is 
over the entire Fermi surface(s), and is weighted by the 
density of states. The zeros mean that the charge life- 
time is infinite; charge is conserved. Linear combinations 
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FIG. 5: Smallest two eigenvalues (absolute values) of the in- 
verse diffuson as a function of momentum. Zeros occur both 
at q = and at q = ±2Q, implying conduction of both the 
spatially uniform spin profile and PSH's. Dashes and dots are 
for the linear and quadratic models seen in Figures ^ and |3j: 
respectively. The Fermi energy Ep has been tuned to match 
the local minima seen in Figure The momentum scale is 
Q = \Q\x, and r Ef = 8. 



of S x , S y have lifetimes 2r/(l ± A /(cos2<9) 2 r + (sin 20) |,); 
spin conduction is obtained only if 

1 = (cos 20)| + (sin 20) 2 F . (3) 

Equation [3] confirms our earlier statement that spin 
conduction can be obtained by tuning for conservation 
of a spin operator cos<j)cr x + sin (j)a y . In this case 
|cos0|, I sin 9 1 are equal to the constants | cos 0| , |sin0| 
and equation [3] is trivially satisfied. Figure [5] shows the 
spectrum of the inverse diffuson in two models which con- 
serve cos (j) a x + sin <p a y . Both models have two null eigen- 
values at q x = 0, corresponding to conduction of charge 
and of S^. Both models also exhibit a single null eigen- 
value at q x = ±2Q, implying PSH conduction. Models 
which both conserve cos <ft <j x +sin <p <j y and have a nesting 
symmetry with 0(Q) = (j) always exhibit PSH's. When 
the nesting symmetry is only approximate, the PSH life- 
time is tpsh = h 2 /^r{(SE) 2 ) F ), SE = g.V k E(Q). The 
gradient in this formula measures violation of the inver- 
sion symmetry Q + k — > Q — k. ^ = k — Q is the displace- 
ment of the Fermi surface from the extremum. 

Equation [3] also confirms that tuning the Fermi energy 
can produce spin conduction. The key is that this equa- 
tion concerns only the Fermi surface: | cos#| and | sin#| 
must be constant there. When Ep is tuned close to a 
local extremum the Fermi surface becomes very small. 
Therefore | cos 9\, \ sin 0\ are nearly constant on the Fermi 
surface, cos (j)cr x + suk/xt^, <p = 6{Q) is nearly conserved 
there, and S<j> conducts freely. 

Tuning the Fermi energy also can produce PSH's. T 
symmetry requires that extrema always come in pairs 
at k = ±Q. The Fermi surfaces S± accompanying 
these pairs possess the nesting symmetry which produces 
PSH's. However there must be no scattering between the 
pair of Fermi surfaces S± and any other Fermi surfaces, 
because cos<pa x + sin^xj^, <j) = 0(Q) will not be con- 



served on the other surfaces. If there are additional Fermi 
surfaces then the disorder potential must be smooth, 
without short- wavelength variations. In this case there 
will be one pair of persistent spin helices for each pair of 
extrema. 

We have computed the spin lifetime when there are 
local extrema. For the spatially uniform spin profile it is 
t s /t = l/((S0) 2 ) Fj SO = <f • V k 0(Q). The PSH lifetime is 
double this value. Our PSH calculation is valid only in 
the diffusive regime where the PSH characteristic length 
I psh = h/2\Q\ is large compared to the scattering length 
l T . SO measures the amount that varies on the Fermi 
surface, because when is constant is conserved and 
r s = 00. £ measures the width of the Fermi surface; when 
Ep approaches the extremum it goes to zero and the spin 
lifetime diverges. For instance, in the quadratic model 
H = + v xk x cr y + Vyk y a x , v y > v x shown in Figure [i] 
the lifetime is r s = r ~^jEj^ F ~y • It diverges when spin is 
conserved (v x = 0) and also when the Fermi energy Ef is 
tuned to the extremum Eq. When the model is tuned for 
rotational symmetry v x —> v y the lifetime becomes very 
small because the local minimum becomes very shallow, 
the Fermi surface stretches along the x axis, and (£ 2 )f oc 
(E F — Eq)/(v 2 — v 2 ) becomes very large. 

Local extrema have already been realized in a 
Bi\- x Sb x , which has six fold symmetry. ARPES 
measurements 2 ^ reveal at least three six-fold degenerate 
minima, with momenta at \Q\ w 0.15, 0.8, 1.1A . How- 
ever the symmetry is too high: substitutional disorder 
causes scattering between all six minima. Moreover the 
bulk gap is very small, and the PSH length scale is so 
short that it may lie in the ballistic regime. 

Local extrema will be found whenever there is an 
avoided band crossing in an anisotropic material. In 
Bii- x Sb x a conventional surface band occurs very close 
to the TI band. Repulsion between these two bands 
causes the observed local minima. Avoided band cross- 
ings have also been observed in very thin TI films - 
the TI bands on each of the film's two surfaces couple 
to each other, causing band repulsion and extremsP^^. 
The remaining necessary ingredient for spin conduction 
is anisotropy. In this respect the recent predictions of 10 
to 1 anisotropy 11 in (3 — Ag^Te and 18 to 1 anisotropy in 
metacinnabapEO are very encouraging. 

When spin conducts - for instance when Ep is tuned 
near an extremum - the magnetoresistance will become 
null or even change sign. If only the spatially uniform 
profile conducts then there will be neither weak local- 
ization nor antilocalization (null magnetoresistance). If 
there are PSH's then there will be a complete reversal 
from weak antilocalization to weak localization, from pos- 
itive to negative magnetoresistance. 

Returning to equation [l] we have calculated the dif- 
fuson operator T>ij which controls spin diffusion via 
T^ij'S = 0. Our calculation is general for all H gen 
but considers only long wavelengths; i.e. momenta near 
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q = 0. For brevity we will present here only the result 
when cos (j)a x + sin (j)a y is conserved. The spin compo- 
nent orthogonal to decays with lifetime r s = r, and 
we have already seen that S z = 0. The spin diffusion 
equation for N and Ss is: 



d t -\V x -D-V x 
-\V X -T d t 



• r 

Vx ■ D ■ V, 





"TV" 







f = ^(V k E(k)) F , D = ^(V k E(k) ® V fe £(fc)) F (4) 

The charge-spin coupling V x • T and the diffusion tensor 
D are determined entirely by the energy dispersion. If the 
Fermi surface is an ellipse with height h and width w then 
f = and D = Sr(E F - E Q ) 2 /H 2 (w~ 2 x(g)X + h- 2 y®y). 
Assuming that the PSH length scale Ipsh = h/2\Q\ ^> 
l T is in the diffusive regime, we have derived the PSH 
diffusion equation: 

(dt + 1/tpsh ~ (V, - 2iQ) • D • (V x - 2iQ))S PSH = (5) 



calculated the spin decay times and spin diffusion equa- 
tions and found two ways to tune for a long spin lifetime 
and spin conduction. The first tuning mechanism is well 
known from quantum wells but new to TIs: tuning the 
Hamiltonian to conserve a spin operator. We found a sec- 
ond tuning mechanism: tuning the Fermi energy near a 
local extremum of the energy dispersion. Neither mecha- 
nism is possible unless the TI surface is anisotropic. Both 
mechanisms cause conduction of a spatially uniform spin 
profile. If the Fermi surfaces exhibit an additional nesting 
symmetry then Persistent Spin Helices will also conduct. 
When spin conduction is realized the TPs magnetoresis- 
tance will be either null or negative, unlike an untuned 
TI where the magnetoresistance is positive. TIs which 
combine anisotropy with avoided band crossings will be 
promising candidates for spin conduction and PSH's. 



The term with two (V^ — 2iQ) derivatives implies that 
small deviations from the spin helix relax diffusively. 



IV. CONCLUSION 

In this article we studied spin conduction in a very gen- 
eral model of TI surfaces with non-magnetic disorder. We 



We acknowledge support from the NSF of China 
(Grant No. NSFC 10876042 and No. NSFC 10874158), 
the 973 program of China (Grant No. 2007CB925000 and 
No. 2011CBA00108)), and the WCU (World Class Uni- 
versity) program of POSTECH through R31-2008-000- 
10059-0, Division of Advanced Materials Science. V. E. 
S. acknowledges the hospitality of AMS, POSTECH. 



* Electronic address: vincent@sacksteder.com 

1 L. Fu and C. L. Kane, Physical Review B 76, 045302 

(2007) . 

2 H. Zhang, C.-X. Lu, X.-L. Qi, X. Dai, Z. Fang, and S.-C. 
Zhang, Nature Physics 5, 438 (2009). 

3 M. Z. Hasan and C. L. Kane, Reviews in Modern Physics 
82, 3045 (2010). 

4 C.-X. Liu, X.-L. Qi, H. J. Zhang, X. Dai, Z. Fang, and 
S.-C. Zhang, Physical Review B 82, 045122 (2010). 

5 L. Fu and C. L. Kane, Physical Review Letters 100, 096407 

(2008) . 

6 J. W. Mclver, D. Hsieh, H. Steinberg, P. Jarillo-Herrero, 
and N. Gedik, Nature Nanotechnology 7, 96 (2011). 

7 J. I. Inoue and A. Tanaka, Physical Review Letters 105, 
017401 (2010). 

8 P. Wenk, M. Yamamoto, J.-i. Ohe, T. Ohtsuki, B. Kramer, 
and S. Kettemann, in Quantum Materials, Lateral Semi- 
conductor Nano structures, Hybrid Systems and Nanocrys- 
tals, edited by D. Heitmann (Springer Berlin Heidelberg, 
2010), NanoScience and Technology, pp. 277-302, U RL 
|http : //dx . doi . org/10 . 1007/978-3-642- 10553- 1.1 l] 

9 A. A. Burkov and D. G. Hawthorn, Physics Review Letters 
105, 066802 (2010). 

10 B. A. Bernevig, J. Orenstein, and S.-C. Zhang, Physical 
Review Letters 97, 236601 (2006). 

11 W. Zhang, R. Yu, W. Feng, Y. Yao, H. Weng, X. Dai, and 
Z. Fang, Physical Review Letters 106, 156808 (2011). 



12 F. Virot, R. Hayn, M. Richter, and J. van den Brink, Phys- 
ical Review Letters 106, 236806 (2011). 

13 R. Egger, A. Zazunov, and A. L. Yeyati, Physics Review 
Letters 105, 136403 (2010). 

14 C.-Y. Moon, J. Han, H. Lee, and H. J. Choi, Physical 
Review B 84, 195425 (2011). 

15 J. Schliemann, J. C. Egues, and D. Loss, Physical Review 
Letters 90, 146801 (2003). 

16 J. D. Koralek, C. P. Weber, J. Orenstein, B. A. Bernevig, 
S.-C. Zhang, S. Mack, and D. D. Awshalom, Nature 458, 
610 (2009). 

17 S. Hikami, A. I. Larkin, and Y. Nagaoka, Prog. Theor. 
Phys. Progress Letters 63, 707 (1980). 

18 H. Suzuura and T. Ando, Journal of the Physical Society 
of Japan 75, 024703 (2006). 

19 E. McCann, K. Kechedzhi, V. I. Fal'ko, H. Suzuura, 
T. Ando, and B. L. Altshuler, Physical Review Letters 
97, 146805 (2006). 

20 A. A. Burkov, A. S. Nunez, and A. H. MacDonald, Physical 
Review B 70, 155308 (2004). 

21 D. Hsieh, D. Qian, L. Wray, Y. Xia, Y. S. Hor, R. J. Cava, 
and M. Z. Hasan, Nature 452, 970 (2008). 

22 H.-J. Zhang, C.-X. Liu, X.-L. Qi, X.-Y. Deng, X. Dai, S.-C. 
Zhang, and Z. Fang, Physical Review B 80, 085307 (2009). 

23 H.-Z. Lu, W.-Y. Shan, W. Yao, Q. Niu, and S.-Q. Shen, 
Physical Review B 81, 115407 (2010). 

24 J. Linder, T. Yokoyama, and A. Sudbo, Physical Review 



6 



B 80, 205401 (2009). 

C.-X. Liu, H. J. Zhang, B. Yan, X.-L. Qi, T. Frauenheim, 
X. Dai, Z. Fang, and S.-C. Zhang, Physical Review B 81, 
041307 (2010). 

T. Oguchi and T. Shishidou, Journal of Physics: Con- 



densed Matter 21, 092001 (2009). 

Highly symmetric TI surfaces (C n v, n > 3) require the 
Rashba Hamiltonian. Lower symmetries^ (for instance 
CW) can produce anisotropic Hamiltonians. 



